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Influence of Coulomb correlations on the quantum well intersubband absorption at 

low temperatures 
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We present a theory for the intersubband absorption including electronic ground-state correlations 
in a doped GaAs/Alo.35Gao.6sAs quantum well system. Focusing on the influence of the Coulomb 
interaction among the carriers at low temperatures, we find that the ground-state correlations lead 
to an increased renormalization and spectral broadening of the absorption spectrum. At T — 1 K, 
its full width at half maximum is increased by up to a factor 3. The inclusion of electron-phonon 
scattering strongly reduces the relative impact of the electronic correlations. 

PACS numbers: 73.21.Fg,78.67.De 



As a source of Terahertz and infrared radiation, semi- 
conductor quantum well intersubband (ISB) transitions 
have been widely investigated experimentally.-^- In par- 
ticular, an understanding of the spectral linewidth is 
of crucial importance for the design and control of 
heterostructure^ and laser devices!^ Although the theo- 
retical description of ISB dynamics has been the subject 
of many studiesp"^ current models for the ISB optics 
become insufficient at low temperatures, yielding an in- 
trinsic absorption linewi dth wh ich is much smaller than 
experimentally observed ! 9 * 11 ! 18 ^ For an equilibrium Fermi 
distribution of the electrons in the ground-state, Pauli 
blocking prevents electron-electron scattering at low tem- 
peratures, while the electron-phonon contribution is too 
weak to explain the experimental findings. To approach 
the experimental findings in the low temperature regime, 
some theoretical descriptions focus on disorder correc- 
tions such as impurities^ or interface roughness P 

In this article, we demonstrate that even in the ab- 
sence of disorder, correlations in the electronic ground- 
state of the doped quantum well yield a line broadening 
at low temperatures not observed before: Coulomb cor- 
relations lead to modified distribution functions which 
partially reduce the strong Pauli blocking in the low 
temperature regim e by o pening up new electron-electron 
scattering channels! 20 ^ 21 ! This leads to additional dephas- 
ing which manifests itself in an enhanced broadening of 
the absorption spectrum. Even if the additional inclu- 
sion of electron-phonon scattering strongly reduces the 
influence of the correlations, a thorough investigation of 
ground-state correlations clearly improves our physical 
understanding of the many-body effects in high quality 
quantum well samples. 

Our approach is as follows: (I) after summarizing the 
overall dynamics of a two-dimensional electron gas in an 
ISB system, (II) we discuss the origin and (III) aspects 
of ground-state correlations and finally show their impact 
on the ISB absorption spectrum. 

(I) ISB dynamics. For our investigations, we consider 
an n-doped GaAs/Alo.35Gao.65As quantum well system 



FIG. 1: (Color online) Sketch of the two energetically lowest 
conduction subbands of an n-doped quantum well. The ar- 
rows symbolize the different Coulomb many-body interactions 
among the electrons. Optical excitation creates an electron 
(a* - ) in band 2 and annihilates an electron (a, - ) in band 

v 2,k,s' y l,k,s' 
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where only the kinetics of the two lowest conduction sub- 
bands are of relevance.^ This is a good approximation 
since we are focusing on low temperatures and subband 
gaps which are typically quite large {k^T < 9 meV <C 100 
meV < £ gap ). To approximate the more realistic finite 
potential of the quantum well, we use an effective well 
width for an infinite potential wellP Non-parabolicity ef- 
fects are included in the form of different effective sub- 
band masses.^ A sketch of the considered quantum well 
in-plane band structure is given in Fig. [II 

The Hamiltonian of the investigated system consists 
of the in-plane kinetics of the confined electrons (Hq e \), 
the Coulomb interaction (Hq), and the semiclassical cou- 
pling to an external light field (H cnl ): 
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Here, a - (a -* ) denotes the creation (annihila- 

i,ki,Si i.ki : Si 



tion) operator for an electron in subband i with an 
in-plane wave vector fej, spin Si, and energy e. t, — 
Si + (h 2 kf/2mi) (cf. Fig. 111). We introduce the com- 
pound index {i} = {i, fej,Sj} to simplify the notation. 
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describes the Coulomb-induced transitions of electrons 
in the states {I, m} to the states {i,j}, where VlJ~ 
is the Coulomb matrix element P The Rabi frequency 
fiem(t) = di2 • ~Ei{t)/h (dipole moment (I12) describes 
the interaction between the external field E(t) and the 
electronic system. The absorption coefficient a(u>) of the 
ISB system is calculated via the complex susceptibility 
x(w) as a(iS) ex wlm[x(w)], x( w ) = P{u)/€qE(uj) with 
the macroscopic polarization P(ui) determined by the mi 
croscopic ISB polarizations p. ■ t 

[cf. Eq. @]P 

The calculations of the ISB dynamics for the po- 
larizations p- ■ t are carried out within a density- 
matrix approach using a correlation expansion! 24 * 25 ^ 
Evaluating the system up to second order, the dy- 
namical equations for the polarizations and the elec- 
tronic populations f ik ~ s = (a r a. ? ) as well as the 

second-order correlations a c {ijlm} = ( a \ i} a \j} a {i} a { m }) ™ 

(( a \i} a {m})( a \j} a {l}) ~ ( a \i} a {l})( a \j} a { m })) are coupled 
by the electron-electron interaction. In the following, we 
will neglect the spin index since the quantities of interest 
are spin independent for our system. 

Evaluating the dynamical equations for the electronic 
polarizations between the two lowest subbands p 12 t with 
the Hamiltonian given above yields 
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The energy £ . g combines the subband energy e. g and the 
energetic renormalization due to the Coulomb exchange 
contribution. The renormalized Rabi frequency tt(t) con- 
tains the Rabi frequency of the external light field £l em (t) 
and a renormalization due t o th e Coulomb exciton and 
depolarization contributions. 8 9 These renormalizations 
result from a Hartree-Fock appr oxima tion, which is first 
order in the Coulomb potential! 26 * 27 ! The second-order 
correlation contribution yields the dephasing functional 
T{p i . £, /. ^} caused by Boltzmann scattering between 
the electrons and includes both diagonal and nondiag- 
onal Coulomb scattering contributions which lead to a 
broadening of the absorption spectrumPSI Depending on 
the nonparabolicity of the bandstructure, a strong can- 
cellation between diagonal and nondiagonal terms can 
occurP For further details concerning the Hartree-Fock 
and kinetic scattering contributions, see Ref. ]T. 

Additional higher-order correlations are included via a 



phenomenological damping 7 of the second-order correla- 
tion functions, 28 which softens the strict energy conser- 
vation in T{p. . g, /. ^} typically used in Markovian ap- 
proaches. In particular, 7 enters T{p i - g, f i g} and the 
correlation correction of the ground-state distribution 
6ft (discussed in the next section) in a way that a strong 
cancellation of the influence of 7 on the spectral broad- 
ening occurs.-^ For the purpose of this paper, it is thus 
justified to focus on a fixed value of hj which is chosen 
as 5 meV here. 

Due to Pauli-blocking in the dephasing functionals 
r{/9.. t, /• t}, the broadening of the spectrum becomes 
very narrow for low temperatures due to the sharp Fermi 
edge. This result, which is not in agreement with experi- 
mental findings, leads us to the assumption that some 
fundamental suppositions in the theory must be cor- 
rected. Here, we address additional ground-state cor- 
relations neglected so far: 

(II) Ground state correlations. In most theoretical de- 
scriptions of £. t, Q, and T, the electronic populations /. t 

in equilibrium are taken to be Fermi distributions / - of 

the non-interacting electron gas! 6 * 8 * 29 ! This yields good 
agreement with the experiment in the high temperature 
regimePSMI At low temperatures, the electron-electron 
ground-state correlations are expected to play an impor- 
tant role for the following reason: While for high temper- 
atures the kinetic energy of the electrons clearly dom- 
inates over the Coulomb repulsion, the latter becomes 
more important for low temperatures where the kinetic 
energy and the Coulomb repulsion may be of a similar 
magnitude. 

To extract the correlation of the ground-state beyond 
the usual second-order Born-Markov approximation, we 
consider a deviation 5f. i — I ■ t ~ f J from the equilib- 
rium Fermi distribution /'J (cf. Refs. 12012X1 for a treat- 
ment of ground-state correlations in metals) and include 
first-order memory effects. The evolution of 5f. t is de- 
scribed by the kinetics flo.el of the non-interacting elec- 
tron gas as well as the Coulomb coupling Hq: 
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We evolve Eq. |5| up to second order in the Coulomb 
coupling similar to the derivation of the ISB polariza- 
tion. Restricting the correlation effects to a single sub- 
band (ground state), we neglect the subband index i in 
the following. In this case, the first-order correlation con- 
tributions vanish, and in second order, we obtain linear 
differential equations for o c <- ., ■. with a time dependent 

^ Z J t Til j- 

inhomogeneity Q(t): 



in ^ 6 h = E [rtfo ?> *') ~ ^ £ k')] V^ , (6) 

k' .q 

7i~a 1/2 = (±Ae + ift 7 )<7j /2 ± W\ % <^'\Q{t), (7) 



with the abbreviations W&*&\ =2V^ - V^'-?'^ and 
Ae = £k_g+£g + ff—Eg ~ £ k> where V denotes the screened 
Coulomb matrix element of the lower subband due to the 
modification of the potential by the presence of the other 
electrons.— Integrating Eq. (|7|, and assuming the mem- 
ory of Q(t) to be small, we can expand the inhomogene- 
ity Q(t') in a perturbation series around the local time t, 
yielding the following expression for a c x ,„: 
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The zeroth-order term ex Q(t) yields the typical Boltz- 
mann scattering contributions which vanish for the Fermi 
distribution functions P3 The second term includes mem- 
ory effects in first order (containing the temporal deriva- 
tive of the source) and leads to 
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which is substituted in Eq. ([6]). 

In a first order iteration, we take the electron occu- 
pations occuring in the inhomogeneity Q(t) to be Fermi 
distributions. This leads to the final expression of the 
deviation 8 ft'- 
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with the abbreviation f- 

J k 



1 



/2. Equation (flO) illus- 



trates the interplay of electrons within one subband to 
renormalize the equilibrium Fermi function: Two elec- 
trons with wave vectors k' + q, fc — <f are annihilated, while 
two electrons with k, k' are created, respectively. Again, 
the phenomenological damping 7 represents higher order 
correlations.^ 

Our result for the correlated ground-state distribution 
function /g is plotted in Fig.plfor three different temper- 
atures: T — 1 K (solid lines)! 50 K (dashed lines), and 
100 K (dotted lines). Comparing ft (dark lines) with 

the equilibrium Fermi distribution fi (light lines), we 

observe a slight decrease (increase) of the electron dis- 
tribution for wave vectors below (above) the Fermi edge, 
already known from electron gases in metals. 20 21 This is 
especially pronounced for T — 1 K, where one now finds 
available states below and a finite population above the 
Fermi edge. Looking at the deviation from the Fermi 
function 8f% (plotted in the inset of Fig. pi), one can see 
the sharp edge around the Fermi energy and a renor- 
malization up to 3% at a temperature of T = 1 K. The 
observed features soften for higher temperatures but re- 
main on the same order of magnitude. 
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FIG. 2: (Color online) Electron distribution function of the 
electronic ground-state including (dark lines) and neglecting 
(light lines) Coulomb correlations for a 5 nm quantum well 
with an electronic doping density of ridop = 6.0 x 10 11 cm~ 2 
and various temperatures. Inset: Corresponding deviation 
Sfr — ft — /i from the Fermi functions. 



Even though the total deviation from the ideal Fermi 
function does not change strongly with temperature, we 
expect the influence of the correlations on scattering pro- 
cesses to decrease for rising temperatures: (i) The rela- 
tive importance of the allowed scattering processes de- 
creases with an increasing softening of the Fermi edge 
(see Fig. [2]) and, at the same time, (ii) the deviation 
from the ideal Fermi function 6ft decreases close to the 
Fermi edge. The scattering processes which go with the 
Coulomb coupling V q ~ 1/|<7| are responsible for the 
heigth of (5/g. Since for larger temperatures the elec- 
tron scattering is less common around the Fermi edge, 
6f% also decreases. 

Even if the calculated 6ft is quantitatively quite 
small, it opens up a new physical scenario for the low- 
temperature regime: 6ft 7^ leads to an occupation 
above and to available states below the Fermi edge, in 
particular for temperatures near K, and allows scatter- 
ing which is otherwise prohibited due to Pauli blocking. 

Besides the dependence on the temperature, the ob- 
served correlation effects also depend on the quantum 
well width, where a slightly stronger deviation is found 
for smaller well widths (not shown). This is due to the 
stronger Coulomb coupling between the electrons for de- 
creasing well widths. The doping density also influences 
the magnitude of 6f%, where an optimal value of the den- 
sity leads to a maximal deviation (cf. also the absorp- 
tion spectra in Fig. [3J). The decrease of the deviation for 
smaller densities can be explained by the short-range na- 
ture of the electronic correlations: for a larger mean-free 
path (Wigner-Seitz radius) r s between the electrons,^ 
the probability of a momentum transfer decreases. On 
the other hand, we expect a decreasing deviation at a 
certain point for large densities since the the kinetic en- 
ergy, which goes with l/r s 2 , overweigths the Coulomb 
interaction, going with l/r s between the electronsPS 
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FIG. 3: (Color online) ISB absorption spectrum of a 5nm 
quantum well at a temperature T = 1 K and different doping 
densities na op (in cm -2 ) including (dark lines) and neglecting 
(light lines) the electronic ground-state correlations. 
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FIG. 4: (Color online) ISB absorption spectrum of a 5 nm 
quantum well with na op — 6.0 x 10 11 cm" 2 including (dark 
lines) and neglecting (light lines) the electronic ground-state 
correlations for different temperatures. Inset: Enlarged view 
of the absorption line shapes. 



(Ill) Absorption spectra. Next, to discuss experimen- 
tal observables, the influence of the ground-state corre- 
lations on the linear quantum well intersubband absorp- 
tion spectrum is studied. Figure [3] shows the spectrum 
calculated from the dynamics of the polarization given 
by Eq. Q for different doping densities at a tempera- 
ture of T = 1 K. Here, the influence of the electronic 
ground-state correlations is clearly visible. We find a 
strong broadening of the absorption line shape includ- 
ing ground-state correlations (dark lines) compared to 
the spectrum neglecting the electronic correlations (light 
lines). The full width at half maximum in the former 
is up to three times larger than in the latter case. In 
addition, the calculation including the correlated ground 
state shows a small energetic renormalization. This sub- 
stantiates our prior claim that even small ground-state 
correlations have a strong influence at low temperatures: 
The absorption linewidth increases due to the generation 



of available scattering states, reflecting a partial cancel- 
lation of the strong Pauli blocking by the electronic cor- 
relations. For the parameters used here, the difference 
between the absorption spectrum including and neglect- 
ing ground-state correlations is maximal for a doping 
density of nd op = 6.0 x 10 11 cm -2 . The shift to lower 
energies for increasing densities is mainly due to the en- 
ergetic renormalization resulting from the Hartree-Fock 
contributions}^ leading to a smaller effective band gap 
between the lower and upper subband. Since we want to 
focus on the impact of the ground-state correlations on 
the absorption spectrum, we will restrict the following 
investigations to a 5 nm quantum well using a doping 
density of nd op = 6.0 x 10 n cm~ 2 . 

In Fig. |4j we show the absorption spectra for differ- 
ent temperatures. Comparing the results including (dark 
lines) and neglecting (light lines) ground-state correla- 
tions, we find that for temperatures higher than 50 K, 
the correlation effects are of no significant relevance as 
discussed earlier in the temperature dependence of 6f% 
(cf . the inset of Fig. [4]) . Instead, nonparabolicity effects 
gain importance, leading to an asymmetric line shape es- 
pecially for smaller well widths due to the occupation of 
higher energetic statesPFor 50 K and below, the absorp- 
tion spectrum shows a significant broadening leading to a 
strongly reduced peak absorption. Furthermore, the cal- 
culation with the correlated ground-state shows a small 
energetic renormalization. 

Although we find that the ground-state correlations 
lead to a significant broadening of the spectrum at low 
temperatures, the linewidth is still strongl y und eresti- 
mated compared to experimental findings! 11 ! 18 ! In or- 
der to present a full and consistent description of ISB 
absorption including all important scattering processes 
which contribute to the absorption line shape, electron- 
phonon interaction must also be taken into accountP^ 
This is assumed to be the major broadening mecha- 
nism for temperatures higher than T =100 KP For 
low temperatures, spontaneous phonon emission yields a 
temperature-independent contribution to the linewidth. 
We therefore extract temperature-dependent Markovian 
electron-phonon scattering rates from earlier microscopic 
calculations and include them as an additional dephasing 
contribution 7 P hon in the calculations™ The resulting cal- 
culated ISB absorption spectra are found in Fig. [5] We 
find that the influence of the ground-state correlations is 
strongly masked by the phonon-induced dephasing. Still, 
a deviation between the absorption spectra neglecting 
ground-state correlations (light solid line) and including 
ground-state correlations (dark solid line) can be found 
for temperatures lower than 50 K, where the the spectra 
neglecting ground-state correlations show less broadening 
than the ones including ground-state correlations. For T 
= 100 K, the difference is hardly visible. 

In conclusion, we have presented a microscopic the- 
ory for the description of ISB quantum well absorption 
including electron-electron contributions and electronic 
ground-state correlations. We showed that by includ- 
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FIG. 5: (Color online) ISB absorption spectrum of a 5 nm 
quantum well including (dark lines) and neglecting (light 
lines) the electronic ground-state correlations, with micro- 
scopically calculated Markovian electron-phonon scattering 
rates for a doping density ridop = 6.0 x 10 11 cm -2 for different 
temperatures. 



ing ground-state correlations, the absorption linewidth 
for temperatures T < 50 K shows significant broadening, 
where the full width at half maximum is increased by 
up to a factor 3 at T = 1 K. The additional inclusion 
of electron-phonon scattering masks the impact of the 
ground-state correlations on the spectral width. In a next 
step, a more consistent treatment of both the electron- 
phonon and electron-electron dephasing, including full 
non-Markovian effectsp^l should be carried out to allow 
deeper insights into the influence of electronic ground- 
state correlations. 
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tric constant £b g = 10.9, longitudinal optical phonon en- 
ergy huiLO = 36 meV, subband gap e gap = 210.66 meV, 
effective masses ml = 0.078mo, m^ = 0.131mo. 
While we include all diagonal terms (proportional to p. . r) 
in the kinetic scattering contributions, we only consider 
the nondiagonal terms in Ref. [9] which are proportional to 
p-tt* since these are the dominant terms counteracting 
the broadening of the diagonal terms. 

The deviation 5/g increases for decreasing 7 while the 
linewidth of the absorption spectrum neglecting correla- 
tion contributions, determined by F{p. . g, f i g}, is reduced 
for decreasing 7. Thus, the broadening effect of 5ft on the 
absorption spectrum is always in a reasonable proportion 
to the reduction of the linewidth so that the influence of 7 



on the spectral broadening is strongly canceled. but finite value for finite values of 7 which is neglected here. 

37 While this term vanishes exactly for 7 = 0, it yields a small 



